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1-S^^ Abstract 

. We perform a comprehensive analysis of the polarization of the top quarks originating 

J> from sbottom-pair production at the LHC, followed by sbottom decays to top+chargino. 

r We study moreover the expected net polarization of top quarks produced in sbottom-to- 

<^ chargino and stop-to-neutralino decays in scenarios with small chargino-neutralino mass 

CNJ difference, where these decays may be hard to distinguish. We show that, in contrast to 

■^ top quarks produced via Standard Model processes, the average polarization of top quarks 

^^ originating from these SUSY decays can obtain any value between +1 and —1. We further 

. study the effect of this polarization on the top quark decay kinematics. On the one hand 

1^ this may be used to construct measures of this polarization, on the other hand it may be 

• ^H used to enhance the search reach in certain scenarios. Exploiting top polarization may 

/\ also prove useful for searches for "natural" SUSY with light higgsinos, which is typically 

j^ very difficult to detect at the LHC. 



1 Introduction 

The discovery of the Higgs boson by both the ATLAS and CMS collaborations [1,2] has provided 
the last missing piece of the Standard Model (SM). Nevertheless there are still fundamental 
problems open in the SM, such as the nature of dark matter, the origin of CP violation, or 
the stability of the electroweak scale, motivating the need for new physics beyond the SM. 
Supersymmetry (SUSY) is one of the best-motivated extensions of the SM for addressing some 
of these issues, and the search for supersymmetric particles is thus one of the primary objectives 
of LHC experiments. 

While SUSY searches at the LHC with ^/s = 7-8 TeV have pushed the mass limits for 
gluinos and light-flavor squarks beyond 1 TeV, current limits on third generation squarks are 
much weaker for several reasons. First, the cross section for direct pair production of stops 
and sbottoms is dominated by gluon-initiated processes and thus drops rapidly as their mass 
increases. Second, this cross section is much smaller than the total squark-gluino cross section 
that contributes to the production of SUSY partners of light quarks. Finally, the decays of 
stops and sbottoms involve some final states with top quarks, which means that the results 
from the generic i^^'^'^+jets searches are not applicable. 

Scenarios with an inverted mass hierarchy, i.e. light third but heavy lst/2nd generation 
squarks, have hence become a new focus of phenomenological studies. In particular it is inter- 
esting to investigate specific signatures of stops and sbottoms as well as methods to determine 
their properties at the LHC. 

Results on stop and sbottom searches in direct pair production are usually quoted in terms 
of Simplified Model Spectra (SMS), which assume 100% branching ratio in a given channel. 
The best limit on stops are currently obtained in the i — ?■ txi channel: based on £ = 21 fb~^ 
of data at 8 TeV, the ATLAS collaboration has excluded stops in the mass range 320-680 GeV 
in the limit of massless neutralinos [3]. Slightly weaker upper limits (600 GeV) are obtained 
from the channel t — > bxt assuming a small mass difference between the chargino and the 
neutralino LSP [4], while sbottom masses below 620 GeV are excluded in the channel b — )• bxi for 
neutralino masses below 150 GeV [5] (both results were obtained by the ATLAS collaboration 
with C = 12.8 fb~^). These limits become much weaker with smaller mass splittings between 
the squark and the neutralino/chargino. Present data therefore allow both top and bottom 
squarks well below the TeV scale. 

For final states involving top quarks, the top polarization can be a useful tool to probe new 
physics at colliders, as it is sensitive to the helicity structure of the production process (for a 
recent summary see e.g. Ref. [6]). There are different ways to measure the top polarization. 
In particular, there is a strong correlation between the polarization of the top quark and the 
angular distributions of its decay leptons. It was shown that this correlation is not affected by 
higher-order corrections [7-9] or new physics contributions to the decay [10-14]. Measures of 
top polarization using angular variables have been constructed [13, 14] The energy fraction of 
the decay leptons may also be exploited as a measure of the top polarization [15]. Moreover, 
it was suggested to make use of boosted top and jet substructure methods in hadronic decays 
to determine the top polarization [16]. Finally, the determination of the top polarization at 
the LHC in ti events using the angular distributions of the top decay products in the top rest 



frame was investigated in [17]. 

In the SM, pair-produced tops are unpolarized, while singly produced tops have polarization 
— 1. In SUSY, the polarization of top quarks produced in the decays of stops or sbottoms can 
take any value between ±1, and may provide information on the underlying SUSY scenario. 
In the context of e^e~ colliders, the top polarization may be used to probe the mixing in the 
stop sector [18, 19]. In this case, one can also use the production cross section which depends 
on the squark mixing. At the LHC, the production cross section is independent of the mixing 
angle, leaving mainly the polarization to extract information on the stop/sbottom mixing (in 
addition to branching ratios, if they can be extracted). 

Possibilities of measuring the top polarization and hence getting information on the stop 
mixing in direct stop-pair production at the LHC with ^/s = 14 TeV were investigated in [20,21]. 
More generally, in [22] some of us showed how the longitudinal polarization of the top quark 
from stop decays into neutralinos, ti — )• txi2) depend on the mixing in both the stop and the 
neutralino sectors, as well as on the mass difference between the stop and the neutralino. More 
precisely, for large mass difference Am = m^ — m^o — rrit, a right-handed (RH) stop produces 
a negative top polarization when it decays into a higgsino and a positive polarization when it 
decays into a bino, and vice-versa for a left-handed (LH) stop. On the other hand, the top 
polarization vanishes in the limit of small Am. Furthermore, as pointed out in [22,23], it will 
affect observables such as the energy distribution of the lepton resulting from the top decay, 
thus impacting the reach for stop searches at the LHC. In fact, different assumptions about the 
top polarization are at least in part responsible for the higher reach in the stop mass exclusion 
by ATLAS (which assumes RH top quarks) than by CMS (which assumes unpolarized top 
quarks) [24,25]. Similar effects are expected for sbottom searches in the top+chargino decay 
channel. 

In this paper, we extend the study of [22] by investigating the behavior of the top polarization 
in sbottom decays into charginos. 

As in the case of stop decays, this polarization may give information on the nature of both 
the sbottom and the chargino. We will show that for a wino-like chargino the polarization 
of the top is always —1, while for a higgsino-like chargino the polarization varies from —1 for 
a RH sbottom to +1 to a LH sbottom. Thus a positive polarization would give indication 
on the higgsino nature of the chargino. Recall that higgsino-like charginos and neutralinos 
have small direct production at the LHC and are therefore hard to probe directly especially 
since the decay products are soft because of the small mass difference between the higgsino 
states. In reality, the situation is more complicated: when the lighter sbottom bi has a large 
LH component, its mass is similar to that of the lighter stop ti (as both are determined by 
the same soft mass parameter), and therefore tops can come from several different channels, 
b — ^ tx~^, t — )■ tx^. At the end, only the net polarization resulting from all different decays will 
be measurable. Relating the top polarization to the underlying properties of stops, sbottoms 
and neutralinos/charginos therefore becomes more challenging. 

The paper is organized as follows. We summarize our notation and conventions in Section 2 
before discussing the fermion polarization in the squark rest frame in Section 3. Polarization- 
dependent observables are presented in Section 4, and the treatment of the boost in Section 5. 



Section 6 contains all numerical results, including the dependence of the polarization on the 
fundamental SUSY parameters, some benchmarks scenarios corresponding to different net po- 
larizations, as well as an analysis of polarization-dependent observables that are relevant for 
searches at the 14 TeV LHC. Section 7 contains our conclusions. 



2 Notation and conventions 

For completeness we review in this section our notation and conventions, and give the relevant 
expressions for sfermion decays into charginos and neutralinos. Overall, we follow the notation 
of [19], which was also used in [26]. 

2.1 Sfermion system 

Ignoring intergenerational mixing, the sfermion mass matrices can be written as a series of 2x2 
matrices, each of which describes sfermions of a specific flavour: 
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with 

m| = M'j^ + m| cos 2/3 (Jg^ — e/ sin^ 9w) + ^/, 

m| = M? + ejm|cos2/3 sin^ ^14/ + rrij, 
a/ = A/ — yU {cot /3, tan/3} , (2) 

for {up, down}-type sfermions; mj, e/ and I^ are the mass, electric charge and the third 
component of the weak isospin of the partner fermion, respectively; M^, M^ and Af are soft 
SUSY-breaking parameters for each family, and /i is the higgsino mass parameter. 

According to eq. (1), A^| is diagonalized by a unitary rotation matrix R^ . The weak eigenstates 
/^ and /^ are thus related to their mass eigenstates /^ and /2 by 

^}]=rHI'\^ R^={ '"'^^, '^^0- (3) 

/s/ V/J \-sindj cos 9 J J ^ ' 

Since the off-diagonal element of A1| is proportional to rrif, this mixing is mostly relevant to 

the third generation, f = t,b, on which we concentrate in the following. The mass eigenvalues 
are given by 

^L = I {m% + ml T ^{ml-miy + Aiajmjy) . (4) 



By convention, we choose fi to be the hghter mass eigenstate, m? < in^^. Notice also that 
m? < rrir < m^. For the mixing angle 67 we choose 

2 2 

cos Of = ""^""^ , sin Oj = "^f' ~ "^^^ . (5) 

^(mj^ - mp2 + (a/m/)2 V ("^/^ ' ^}J^ + («/"^/)^ 

The /l^/r mixing is large if (m| — m^ ) < {afirif), with |cos^?| > 4= if m? < m? and 
I cos^?| < 4= if rrif < m? . 

2.2 Neutralino system 

In the basis 

^; = (-<-^A^<,v^°,J (6) 

the neutralino mass matrix is: 

/ Ml —iTiz sin 6'vi/ cos (3 mz sin ^vi^ sin /? \ 

. , _ M2 mz cos ^VK cos /3 — m^ cos ^vi^ sin /3 , . 

—mz sin 6w cos l3 mz cos 6w cos l3 —fj, 

\ m^z sin 9w sin /3 — m^ cos 9w sin /3 — /x J 

The matrix of eq. (7) is diagonalized by the unitary mixing matrix A^: 

NM.jsiN'^ = diag(m^o, m^o, m^o, m^o) , (8) 

where m^o, n = 1,...,4, are the (non-negative) masses of the physical neutralino states with 
m^o < .... < m^o. 

2.3 Chargino system 

The chargino mass matrix is: 

\ y/Im^wcosp n J 

It is diagonalized by two unitary matrices U and V, 

UMcV^ = diag(m^±, m^±) , (10) 

where m-± are the masses of the physical chargino states with m,-± < m-±. 



2.4 Sfermion interaction with charginos and neutralinos 

The sfermion interactions with charginos, which will define the b — )■ tx'^ decays, are 

£^,;^± = gu{lf^Pj, + kf^PL)xtd. + 9d{ltjPR + ki;PL)xt'u, + h.c. (11) 

where i,j = 1,2, u (u) stands for up-type (s)quark, and d (d) stands for down-type (s)quark, 
g is the SU(2) coupling constant. The coupling matrices l-^ and A;-^ are 

llj = -VjiR\^ + ht V,-2-R*2 1 hj = -UjiRii + h Uj2Ri2 , (12) 

jJ^ = h,U,2li,, kl = h,V,2R\,^ (13) 

for stops and sbottoms, with the Yukawa couplings hf given by 

u ^* h ^^ (^A\ 

^* = ~7^ ^~R ' ^' = ~7^ ~R ■ ^^^^ 

VZmwsmp \'2mwCosp 

The sfermion interactions with neutralinos are (? = 1, 2; n = 1, ...,4) 

%x» = 9fifLPR + hLPL)xUL + 9fihLPR + fLPL)x^fR + ^-C- 

= 9f{aLPR + hLP^x^l + \i.c. (15) 



(16) 
(17) 



fL= ^2^Nn2 - \ta.newNni) , 
/i = -^tan^H'iV„i, 
4„ = -/i,iV„3 = 4; (18) 

for stops (/ = t) on the left and sbottoms (/ = h) on the right. 

3 Fermion polarization in the sfermion rest frame 

For a sfermion decaying into a chargino or a neutralino, the gaugino interaction conserves the 
helicity of the sfermion while the higgsino interaction fiips it. We define the average polarisation 
of the produced fermions as 

where o-(±, ±) is the cross section for a positive or negative helicity fermion respectively. In 
general one expects non trivial polarization effects only for the third generation fermions, fur- 
thermore only the polarization of the top quark or the r can be measured. In this paper we 
will discuss only the top polarization. 
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= -^{Nn2 + \ta.new 
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= ^tan9wN^i, 




h^Rn 


= -htNnA = h'^n, 





3.1 h —^ tx decays 

For tops coming from bi ^ txj decays, the polarisation given by eq. (19), is 

- 2 - 2 

^* = , ^f'^^ ~^'^^ ^/^, , (20) 

V^ij) + ('ij) ~ 2 Kjj l,ij J 2 

with the factors /i and /2, purely kinematical in origin, given by 

{P^- ■ St) m^- 

/i = rut—^ r, /2 = mt- — r. (21) 

[Pt ■ P^j) {Pt ■ P^j) 

Here mt, pt and St denote the top mass, momentum and longitudinal spin vector, respectively, 
and p^- is the momentum of the chargino. In the rest frame of the decaying sbottom, these 
factors become 

Ai(m?,m2,m2 ) 2mtm-- 

/l ~ 2 2 2 ' ^2 — — 2 2 2 ' l^^/ 

TTij — m _ — 1717 mi — mt- — m,f 

bi Xj h Xj 

with the function A defined as 

\{x,y,z) = x'^ + y^ + z'^ — 2xy — 2yz — 2xz. (23) 

The factors /i — ?> 1 and /2 — ?■ if mt was negligible. Moreover, for large Am = mr, —m,-- —mt, 

_ I Xj 

/i — 7> 1 irrespective of mt. For hi — )■ tx' decays, we have 

K)-{l\,) = h'tV^2Cos^di-ihU,2sm9-,-U,iCos9if 

= {hlVf^ - U^) cos' % - hlUf^ sin' 9-, + h U,iU,2 sin 29~, . (24) 

r 2 r 2 

For ^2 decays, the corresponding expression (kL) — (/2,) is given by the RHS of eq. (24) with 
cos' 9^, sin' 9f^ interchanged, and a change in sign of the term oc sin 2^^. 

It is interesting to consider certain limiting cases for the bi — > txi decay. For M2 ^ \n\, 
the chargino is wino-like, and the mixing matrices U and V are given by:^ 

The resulting top polarization Vt is then: 

Vt = ^%^ = -h . (26) 



^Up to an overall phase. 



Therefore, for large enough Am, we expect P^ ^ — 1 in case of a pure wino-hke chargino, ir- 
respective of the sbottom mixing angle. (It should be noted however that in practice for the 
range of parameters we will be considering all entries of the matrices U and V are non-zero and 
will be taken into account.) 

For a higgsino-like chargino, i.e. \fi\ ^ M2, the mixing matrices U and V approach 
and the resulting top polarization is 

* [{ht COS e-,y + {h sin e-,y-hthsm2e-j2]' ^ ' 

Therefore, for a higgsino-like chargino, the top polarization depends on the sbottom mixing. 
In the limit of pure LH or RH sbottoms we have 

II : cos^s = l, Vt^+fi, (29) 

bn : cos9~, = 0, V,^-fi. (30) 

In other words, a wino-like chargino, whose interaction conserves chirality, couples to a left- 
chiral sbottom and a left-handed top (recall that just like the W^, a wino has only left-chiral 
interactions) thus always giving a top polarization close to —1, similar to single top production 
in the SM. The higgsino interaction, on the other hand, flips chirality and couples a right (left) 
chiral sbottom to a left (right) handed top; the top polarization thus can vary from —1 to 
-|-1 depending on the sbottom mixing angle and the bottom Yukawa coupling, which becomes 
relevant at large tan/3. This will be discussed in more detail in Section 6. 

3.2 i ^tx^ decays 

For completeness we also summarize the case of stop decays into neutralinos, c.f. [26]. Analogous 
to eq. (20), the average top polarization from tj — )■ tx° is given by 

(bin) +(o4) -2 6LaL/2 
with the obvious replacement p-- — )■ p^o in eq. (21). Using eqs. (16), (17) and (3), the stop 
couplings are written as 

7=2 i""-^ I 



a\n = p Nn2 + - tan 9wNni COS Of - htNni sin 9i 



2\f2 
6*i„ = -—iwiQw^nx^-y^Qi-htNl^cosQi (32) 

o 

Substituting these in the expression for the polarization, eq. (31), one easily sees that in the 

case of large mass differences between the stop and the neutralino (/i — ?> 1) a RH stop will lead 

to P^ = — 1 when it decays into a higgsino and to P^ = +1 when the decay is into a gaugino; a 

LH stop will lead to the opposite polarizations. 



4 Effect of top polarization on the decay kinematics 

The V — A interaction involved in the decay of the top quark, t — )■ Wb -^ ii' b, where i and i' 
denote the decay products of the W, imphes definite correlations between the direction of top 
spin and the top decay products. These are most clearly understood in the rest frame of the 
top quark. Since the top decays before it hadronises, these correlations are not washed out by 
the hadronisation process. We discuss explicitly only the case of the top, given that top and 
anti-top can be distinguished by the charge of the decay lepton. Note also that we neglect the 
effect of the off-diagonal elements of the CKM matrix and use BR{t — )■ Wb) = 1. 

Consider a top quark ensemble with degree of polarization Vt- In the top quark rest frame, 
the angular distribution of the decay product / is given by: 

[1 + KfVt cos 9 forest) , (33) 



Ttdcos9fre.t 2 



/,rest 

where F^ is the partial decay width of the top and 6 forest denotes the angle between the mo- 
mentum of the / and the top spin vector. Hf is called analyzing power of the decay product 
/ [27]. It is 1 for a positively charged lepton or a li quark. For a u quark or ui, Ku,u = —0.31, 
while for b and W the values are Kb = —0.4 and kw = 0.4, respectively. The maximal value of 
Ki = 1 means that the charged lepton is the most efficient polarisation analyzer. Corrections to 
the values of Kf can originate from any deviation of the tbW coupling from the standard V — A 
structure and/or from higher order QCD and QED corrections. The leading QCD corrections 
to Kb, Kd and Ku are of the order of a few percent, somewhat decreasing its magnitude [9]. The 
values of ki and Kd on the other hand do not receive any corrections from the anomalous tbW 
coupling at leading order [14] . Hence the angular distribution of the decay lepton in the rest 
frame reflects the polarization of the decaying top quark faithfully, even in the presence of such 
anomalous couplings. Thus it is an unambiguous measure of top polarization effects. 

It is also interesting to consider how polarization affects the kinematic distributions of 
the top decay products in the laboratory frame. This allows to construct useful polarization- 
dependent observables. The use of laboratory observables to measure top polarization would 
obviate the need for reconstruction of the top rest frame. This is desirable as such a recon- 
struction may not be always possible. The correlation between the polarization of the top and 
the different kinematic variables of the decay products can be obtained from eq. (33) and ap- 
propriate Lorentz transformations. Just like the angular distribution in the top rest frame, the 
energy integrated decay lepton angular distributions in the laboratory frame are also unaltered 
to linear order in the anomalous tbW coupling [10-14]. For all the other distributions, including 
the energy distribution of the decay lepton, a deviation of the distribution from the unpolarized 
case can not be uniquely attributed to the top polarization. 

While constructing polarization dependent observables, it is worth recalling that the decay 
product distributions in the lab frame are influenced not only by the top quark polarization, 
but also by the boost Pt from the top-quark rest frame to the laboratory frame and by the 
transverse momentum p^ of the top quark. Here we will use a boost parameter based on the 



total momentum of the top \'pt\ and the top energy Et 

A = f . (34) 

As an example we consider the lab-frame polar angle di of the lepton w.r.t. the top quark 
direction. Due to the top boost, 9i is smaller than its counterpart in the rest frame ^i,rcst- Thus, 
the distribution of 9i in the lab frame is more strongly peaked towards for a stronger top 
boost as well as for a more positively polarized top quark. One can then define a polar angle 
asymmetry Aq^ as 

_ ajOi < n/A) - ajOi > n/A) 
'' a{6i<'K/^) + a{ei>'K/A)- ^ ' 

In addition to the polar angle, one can study the azimuthal angle distribution. To this 
end, we choose the proton beam direction as the z-axis and define the top production plane 
as the X ~ z plane. Moreover, we identify the positive x component with the direction of the 
top quark. At the LHC, since the initial state has identical particles, the 2;-axis can point in 
the direction of either proton. This symmetry implies that one cannot distinguish between an 
azimuthal angle and an angle 27r — 0. In the rest frame this variable does not depend on the 
longitudinal polarization, but in the lab frame it picks up a dependence on ^;,rest through the 
top boost. For positively polarized tops it is peaked at 0^ = and 0; = 27r, with a minimum 
at 0; = TT [13, 14]. It should be noted that nonzero pj also causes the (pi distributions to peak 
near 0^ = and 0^ = 27r, independent of the polarization state of the t quark. In other words, 
the peaking at 0; = and 2ii is caused by kinematic effects, even for an unpolarized top. It is 
enhanced even further for a positively polarized top. For a completely negatively polarized top, 
the pure polarization dependent effects can sometimes even overcome the peaking caused by 
kinematical effects. The peaks of the distribution then shift a little away from = and 2tt. 
More importantly they lie below those expected for the positively polarized and unpolarized 
top. The relative number of leptons near = and 27r is thus reduced progressively as we 
go from a positively polarized to unpolarized to a negatively polarized top. For normalized 
distributions the ordering is exactly the opposite at = vr where the relative number of leptons 
increases as we go from a positively polarized top to a negatively polarized top. This shape 
then motivates the definition of the azimuthal angle asymmetry [14]: 

_ O-(COS0; > 0) - q-(cOS0; < 0) 

^' ~ a(cos 0/ > 0) + a(cos 0^ < 0) ' ^ ^ 

where a is the fully integrated cross section. Note that a higher top polarization or a stronger 
top boost will result in a more sharply peaked (pi distribution and thus yield a larger asymmetry. 
It is also useful to consider energy observables. Although they are not completely indepen- 
dent of an anomalous thW coupling as mentioned above, they do carry information about the 
top polarization. In fact, for a positively polarized top, the energy and the transverse momen- 
tum distributions for the lepton are shifted to higher values as compared to the unpolarized 
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or negatively polarized case. Since top quarks produced from SM processes are either unpolar- 
ized or negatively polarized, this feature of the Ei and pip distributions for positively polarized 
top quarks can provide an effective discrimination against the SM background. Since the Kb 
and Ki have opposite signs, the effect of top polarization on the energy and px distributions 
of the 6-jet in the laboratory frame is exactly in the opposite direction to that for the lepton 
distributions. In Ref. [28], this feature was employed in constructing a discriminator of top 
quark polarization using the Pt of the 6-quark. Furthermore, the energy distribution can be 
of particular use when the top quarks are highly boosted. In this case, the effect of the boost 
on the angular distribution may mask the polarization and an accurate determination of the 
angles (for asymmetries) may be difficult. It was shown in [29] that in a kinematic regime 
where the tops are heavily boosted, the ratios 

Eb El 

z = — -, u = — — , [67) 

Et El + Eb 

are sensitive to the polarization state of the top quark. Here Et, Eb and Ei are respectively the 
(lab frame) energies of the top quark, and of the b quark and lepton coming from its decay. The 
analysis of [29] was at the LO parton level, but in practical applications one may also consider 
Eb to be the energy of e.g. a b jet. Note that the ranges of z and u are given in principle by 

0<z,u<l, (38) 

although there will be a cut-off at high and low values due to the finite b quark and W boson 
masses. In the coUinear limit (3t = 1, the normalized distribution j^^ can be in fact computed 
analytically [29]. It is peaked at lower values of z for a positively polarized tops, and at high 
values of z for negatively polarized tops. In case of the u distribution, which has to be computed 
numerically, even in the /3i = 1 limit, the peak is shifted by about 0.1 for Vt = —I compared 
to the unpolarized case; whereas for Vt = 1 the normalized distribution is weighted towards 
larger values of u. One can of course define these observables for any value of a cut on the 
top boost. However, at low boost values, both z and u are increasingly contaminated with 
contributions that are independent of Vt, thus reducing their effectiveness as discriminators of 
top polarization and/or new physics parameters. We will show later how these distributions 
may be exploited for quantitative measures of polarization. 
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5 Boost treatment 

Before proceeding to the numerical analysis, we find it useful to give some details on the 
treatment of the boost. Note that the overall boost of the top quark in the laboratory frame 
depends on the boost of the sfermion, that of the top in the sfermion rest frame and the angle 
of emission of the top with respect to the sfermion. As a consequence, the relation between 
the top polarization in the sfermion rest frame, discussed in Sections 3 and 6.1, and the top 
polarization measured in the lab frame depends on all these. 

We explain the procedure we have used to obtain the top quark polarization in the lab frame 
with the concrete example of 6i — )■ txi decays. The same treatment will later be applied to all 
other decays, including stop decays to top+neutralino, in order to obtain the total polarization 
from all sbottom and stop decays into tops, i.e. the quantity that is actually relevant for 
experiment. 

First, we define our frame of reference such that the decaying sbottom is at rest and the 
top momentum lies in the x — z plane. In this frame, the momentum vectors of the sbottom, 
top and chargino are defined as: 

p-,^ = K,, 0,0,0) (39) 

Pt = {Et,\pt\smi^,0,\pt\cosi^) = {Et,pt^,0,PtJ 
Pxi = (^^-,-|pt|sin^,0,-|pt|cos^) 

where, f? is the angle between decaying top and chargino and \pt\ is the magnitude of the top 
momentum given by: 

A5(m? ,m?,m2±) 
\Pt\ = i ^ (40) 



Furthermore, we specify the z-component of the top spin vector as: 



.2 



^3 _ I \Pt\ Et Pt, ^ Et Pt, 



5 = — , -rT,0,^f^ (41) 

\mt rrit \pt\ rrit \pt\J 

The top-polarization in terms of helicity amplitude formalism is defined as: 

f Ml^lrfcos^- r^^^^rfcos^ 

J a cos t) J a cos v 



V, 



J a cos w J a cos t) 

jkfj-llj) I^tip^- ■ sf)d cos^ 

(kfj + Ifj) J{Pt ■ p^-)dcos'd - 2kij lij J mt m^-dcos'd 



(42) 



Now we boost the system in the positive z-direction. Note that the spin vector s^ is not a 
Lorentz vector, hence the dot product {p^- ■ s^) is not Lorentz invariant. Thus, the resulting 
top polarization has a dependence on the boost. 

^Since we are interested in the longitudinal polarization of the top, wc are not concerned about the spin 
vectors in the x and y directions. 
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6 Numerical analysis 

For the numerical analysis, we choose to work within the framework of the general R-parity 
and CP conserving MSSM with parameters defined at the electroweak scale. The relevant soft 
terms for our analysis are the left and right 3rd generation soft masses Mq , Mq , M^ and 
the trilinear couplings At, Af, entering the stop and sbottom mass matrices, together with the 
gaugino masses Mi and M2, the higgsino mass parameter /i, and tan/3. The top quark mass is 
fixed at 173.2 GeV. 

We consider the two cases of higgsino-like or wino-like xf- For the example of a higgsino- 
like chargino, we set /i = 350 GeV and M2 = 1000 GeV; the neutralino sector is fixed by 
Ml = 500 GeV. This gives m^± ~ 352 GeV, m^o ~ 358 GeV and m^o ~ 343 GeV (for 
tan/3 = 10, but showing only little variation with tan/3). For the wino-like case, we reverse 
the parameters, setting M2 = 350 GeV and ^u = 1000 GeV. With Mi = 250 GeV, this gives 
m-± ~ m^o ~ 360 GeV and m^o ^ 247 GeV. ^ 

In general, we use Sof tSUSY [30] to compute the full MSSM spectrum. However, when 
considering only the bi — )■ Xit decay, we directly take the lighter sbottom mass, m^ , and the 
sbottom mixing angle, cos 61, as free parameters; the exact values of the stop/sbottom soft terms 
are not necessary in this case. For computing sparticle decays, checking flavor observables, etc., 
we use micrOMEGAs [31,32]. Cross sections are computed at NLO with Prospino [33]. Finally, 
for generating distributions, we use MadGraph [34, 35] with v4 model files. 

6.1 Parameter dependencies for 61 

Let us first discuss the results for top polarization in the sbottom rest frame. Figure 1 shows 
the polarization of the top quark coming from the 61 — )■ Xit decay as a function of the sbottom 
mixing angle, for the case of large Am; concretely we take m^^ = 685 GeV and, as mentioned 
above, /x = 350 GeV and M2 = 1000 GeV (M2 = 350 GeV and^ = 1000 GeV) for the higgsino 
(wino) case, with Mi = 500 GeV. As discussed in Section 3.1, for a wino-like chargino the 
polarization of the top is always ?^ — 1 regardless of the nature of the decaying sbottom. For 
a higgsino-like chargino, on the other hand, V^ varies from —1 for a bi to +1 for a 6r. For 
small tan/3, i.e. small bottom Yukawa coupling, the transition happens very quickly around 
cos6i ~ 0.1 — 0.25. For large tan/3, where top and bottom Yukawa couplings are equally 
important, the transition happens more slowly and for maximally mixed sbottoms; as a result, 
all values oi V^ G [— 1,+1] can be easily obtained without having to fine-tune the value of 
the sbottom mixing angle. We also note that for cos9i = 0, V^ = —0.92; when cos6'f, slightly 
increases above zero, the top Yukawa coupling starts playing a role in the denominator of 
eq. (28), causing V^ to first decrease to —1 before increasing to +1 for larger cos^^,. A similar 
effect can be seen for cos^j, ~ 1: V^ = +1 is obtained for cos^j, = 0.7, while at cos9i = 1, 
Vt = 0.945. 



^When the rest of the SUSY spectrum is needed, we set the slepton soft terms to M^, = Mj=,, — 500 GeV 
(i = 1...3), the squark soft terms of the first two generations to Ma — Mjy, = Mj^ — 2 TeV (j = 1,2), the 
gluino soft mass M3 = 1.5 TeV, and the pseudoscalar Higgs mass mA = 1.5 TeV. The right stop mass parameter 
Mjj and the trihnear couphng At are adjusted such that mu ~ 126 GeV. 
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Figure 1: Top polarization in 6i — )■ X\t decays as a function of the sbottom mixing angle, for 
nil ~ 68^ GeV. For the higgsino-like case /i = 350 GeV and Mi = 1000 GeV (in red), while for 
the wino-like case M2 = 350 GeV and /i = 1000 GeV (in blue). The solid, dotted, and dashed 
lines are for tan/3 = 10, 30, and 50, respectively. 

In order to understand the dependence of V^ on the relative strength of the Yukawa cou- 
plings, let us simplify eq. (20) for a completely higgsino-like chargino, cf. eq. (28). We get 



V, 



sin-^ 6»r 



"■b 



h 



COS^ Or , h'i 

b I b 

sin^ Or ~'~ h? 



cose^h^ „ 
"^ siner ht •''^ 



(43) 



Now recall that /i ~ 1 away from the kinematic boundary, while fi — ?■ 0. Hence, for 
cos^ ^l^/sin^ Oi > hl/hl, V^ > and for cos^ %/ sin^ 6^ < hl/hf, V^ < 0. This effect is illustrated 
in Fig. 2. Also, for a fixed value of cos 6*^, |Pj| decreases with increasing tan/3. 

The dependence on the available phase space. Am, is illustrated in Fig. 3. The left panel 
shows the case of a higgsino-like chargino while the right panel shows the wino-like case. The 
top polarization as described in eq. (20) is directly proportional to /i, which goes to zero when 
Am — !• and, as a result, V^ — )■ 0. This can very well be seen in both the higgsino-like 
and the wino-like scenarios. Away from the kinematic boundary, /i — )■ 1 and fi — ;■ 0. For 
Am > 100 GeV, /i is very near 1 and the polarization approaches its limiting value. Here, 
the contribution of the fi eq. (28) is somewhat suppressed by the ratio of Yukawa couplings or 
the sbottom mixing angle. Note also that in the higgsino case, for certain values of the mixing 
angle, e.g. cosOi = 0.5, increasing the value of tan/3 can flip the sign of the polarization, this 
corresponds to the transition cos^ 6'^/ sin^ 6'^ = hl/h'f. Contours for the top polarization in the 
cos 61 vs tan/3 plane are displayed in Fig. 4 for both small and large values of Am. In both 
cases any value of Pt can be reached for any value of tan/3, although for Am, = 10 GeV and 
small tan /3 the polarization changes very rapidly with the sbottom mixing angle as a result of 
the kinematic factor fi in the denominator of eq. (28). 

To see whether there is any hope to actually measure this polarization, we also need to 
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Figure 2: Dependence of the top polarization in 61 — t- x^ t decays on the relative strength of 
the t and b Yukawa couplings and the sbottom mixing angle. 

consider the sbottom production cross section and decay branching ratios. The branching 
ratios for 61 decaying into charginos and neutralinos as function of the sbottom mixing angle 
are depicted in Fig. 5. The plot on the left shows the higgsino-like case. At small tan/3, the 
decay into Xit is always dominant even though the branching ratio is reduced when 61 ^ Br 
(cos 6*^ ~ 0), because the coupling strength of a 6^ to a higgsino is proportional to the bottom 
Yukawa coupling. Therefore, at small cos^^;, sbottom decays into charginos and each of the 
three lightest neutralinos become comparable (about 20-30%). For tan/3 = 50, the bottom 
Yukawa coupling is large and the branching ratio to Xit is about 40-60% over the whole cos 6'^ 
range. 

In the right panel of Fig. 5 the Xi ^^id x? are wino-like, X2 is bino-like and X3, X4 are higgsino- 
like. The wino states do not couple to bji, so for cos 6'^ ~ the decay 61 — )■ X2^ dominates. As 
the hi admixture increases, the branching ratio into the wino-like chargino increases because 
the coupling is proportional to top Yukawa coupling. The BR(6i — )■ x'^t) saturates at around 
60% because it has to compete with the bi — )• Xib decay, which has a BR around 40%. 

The cross sections for sbottom (and stop) pair-production as a function of the sbottom 
(stop) mass, calculated at NLO at ^/s = 14 TeV with Prospino [33], are shown in Fig. 6. 
The processes pp — )■ b^b*, and likewise i-i*, {i = 1,2) proceed through pure QCD interaction, 
the mixing angle in the stop or the sbottom sector enters the cross section calculations only 
through 0{as) corrections involving ttg and the four-squark couplings. Thus the cross section 
for stop-pair production is very similar to that of sbottom-pair production. Mixed pairs tit2 
and 6162 cannot be produced at lowest order since the gtt and ggtt vertices are diagonal in 
the chiral as well as in the mass basis, iibi production cross sections are suppressed by at 
least 2 orders of magnitude as compared to tjtj and bibi productions, we ignore this production 
mode for stop and sbottom. The production cross section changes rapidly in the interesting 
mass range between 500 GeV and 1 TeV. For masses around 650 GeV, as considered in the 
benchmark scenarios below, we find a cross section of about 150 fb, which quickly falls off to 
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Figure 3: Top polarization in bi — ?■ Xi ^ decays as a function of Am = nii^ — m-~ — rrit for 
various choices of cos ^5; on the left for fi = 350 GeV and M2 = 1000 GeV (higgsino case), on 
the right for M2 = 350 GeV and /i = 1000 GeV (wino case). Solid lines represent tan/3 = 10 
and dashed lines are for tan/3 = 50. 
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Figure 4: Contours of top polarization in the cos^ versus tan/3 plane for the higgsino case with 
\i = 350 GeV and M2 = 1000 GeV, for a fixed mass difference of Am = 10 GeV (left) and 
Am = 200 GeV (right). 

below 50 fb for around 800 GeV, see the small inset in Fig. 6. 
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Figure 5: Branching ratios of 6i as a function of cos6j^ for tti^; = 685 GeV, on the left for the 
higgsino case with fi = 350 GeV and M2 = 1000 GeV, on the right for the wino case with 
/i = 1000 GeV and M2 = 350 GeV. The sohd (dashed) hues represent tan/3 = 10 (50). The 
remaining parameters are as for benchmark point bm-2 in Table 1. 




Figure 6: Stop and sbottom pair production cross sections at the LHC with -y/s = 14 TeV, 
computed at NLO with Prospino. 
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6.2 Benchmark scenarios 

In a complete MSSM scenario, tops can also come from decays of stops, ti^2, and of the heavier 
sbottom, 62- Some of these decays can be distinguished by their different signatures, for instance 
^1 ~^ ^ xr ~^ ^^ Xi versus t^ ^ t X2 -^ t ^X? or t^ — )■ t Xi- Others, like b-^^ and 62 both decaying 
to t Xi give identical signatures. It is hence necessary to consistently add up the polarizations 
resulting from all processes which cannot be distinguished in the experimental analysis. In 
particular, if the masses of the xf^ X2 ^^^ Xi ^^^ close to each other, as is typically the case 
in higgsino or wino LSP scenarios, the xf and X2 decays into the LSP will lead to soft decay 
products which are likely to be missed. In such a case, the processes 

pp -^ bj>l, W^txi 
pp -^ h^hl, h^^txi 
pp ^ tiil k^txh (44) 

may all contribute to the ti + EJ^^^^ signature, and the total or "net" top polarization relevant 
for the experimental analysis will be a result of the stop and sbottom production cross sections, 
their decay branching ratios and the relevant boosts for going from the rest frame to the lab 
frame. (We neglect the ^2 in eq. (44) because, in order to achieve rrih ~ 126 GeV, at least one 
of the stops should be heavy and will thus have a very low production cross section.) 

In this context, we remind the reader that the situation of only bi being light and all other 
3rd generation squarks heavy only occurs for bi ~ Sr. For 61 ~ bL, also the ti will be close in 
mass (or even lighter because of L-R mixing), because both their masses are governed by the 
same mass parameter Mq . 

For illustration and to allow a complete analysis, we present in Table 1 a set of 7 benchmark 
points which exemplify different scenarios of stop/sbottom mass and mixing patterns and re- 
sulting top polarizations. The production cross sections, branching ratios and top polarizations 
Vt originating from different decay processes in the respective squark rest frame are listed in 
Table 2. Moreover, Table 2 gives the net polarization Vt in the laboratory frame, together with 
the values for the polar angle and azimuthal angle asymmetries defined in eqs. (35) and (36), 
summing over all processes that cannot be distinguished from b^ — )■ tx]".^ 

Since the dependence of the top polarization on the sbottom and stop masses and mixings is 
most interesting for the higgsino scenario, most of our examples focus on this case. Concretely, 
points bm-1 to bm-4 and bm-6 have /i = 350 GeV and M2 = 2Mi = 1000 GeV, as used earlier 
in this paper, leading to a 97% higgsino LSP and a x^^ — Xi mass difference of about 9 GeV; 
bm-5 has a similar electroweak-ino pattern but for a lighter mass scale. The case of Xi ~ W^^ is 
exemplified in bm-7. For all points, the parameters in the squark sector are adjusted such that 
mi ~ 650 GeV (to avoid kinematic effects on Vt) and m^ ^ 126 GeV. All our benchmark points 
are for tan/3 = 10; we do not present specific points for large tan/3 because they would not add 
any new features w.r.t. points bm-1 to bin-7. Note also that we have chosen our benchmark 
points such that they lie (just) outside the current exclusion limits — prospects should be good 



^Of course dedicated simulations would be necessary to decide whether or not specific processes can be 
distinguished. This is however beyond the scope of this paper. Instead we generally treat 61 — > t^T ^^"^ 
ii ^ t Xi 2 ^^ indistinguishable if the xf — Xi mass difference is small, below about 20 GeV. 
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bm-1 


bm-2 


bm-3 


bin-4 


bm-5 


bni-6 


bm-7 




571 
1500 
2200 


1300 

572 

1300 


630 

597 
2200 


610 

599 

2200 


582 
1500 
2200 


850 

601.5 

1500 


600 

609 

2200 


Ml 
M2 


500 

1000 

350 


500 

1000 

350 


500 

1000 

350 


500 

1000 

350 


200 
400 
150 


500 

1000 

350 


250 

350 

1000 


cos 61 


650.41 

1537.39 

0.999 


650.44 

1329.77 

0.006 


650.85 

707.36 

0.10 


650.13 

688.72 

0.15 


650.08 

1538.10 

0.999 


650.17 

884.47 

0.02 


650.71 

692.45 

0.42 


cos Of 


630.18 

2210.95 

0.995 


1236.20 

1560.58 

0.85 


687.92 

2208.66 

0.995 


668.37 

2209.42 

0.995 


633.84 

2212.19 

0.996 


818.99 

1518.19 

0.96 


667.61 

2205.97 

0.996 


^xr 

^x? 
^x^ 
^x§ 
^x2 


351.91 
1013.41 
343.06 
357.41 
502.99 
1013.95 


352.73 
1021.05 
343.39 
358.22 
500.89 
1021.33 


352.04 
1013.45 
343.10 
357.58 
502.32 
1014.02 


351.98 
1013.42 
343.08 
357.50 
502.45 
1013.98 


144.16 
430.32 
126.52 
159.18 
213.20 
431.95 


352.05 
1016.09 
343.13 
357.22 
502.61 
1016.63 


360.35 

1008.07 

247.08 

360.10 

1002.10 

1008.24 




3.53 
3.04 


3.61 
3.04 


3.51 
3.03 


3.51 
3.04 


3.26 
3.04 


3.78 
3.03 


3.31 
3.03 



Table 1: Parameters and masses (in GeV) for seven illustrative benchmarch points. All points 
have tan/3 = 10 and A^ = 100 GeV, while At ^ 3 TeV is adjusted such that rrih ~ 126 GeV. 
The values for BR(i? — )■ Xg'-f) and BK{Bs — )■ fi^fj^~) are given in units of 10^^ and 10~^, 
respectively. 

to test them early in the next phase of LHC running at 13-14 TeV. The detailed characteristics 
of the various benchmark points are as follows: 

• bm-1 features a higgsino scenario with an almost pure LH bi, for which we expect a top 
polarization in &i — )■ txr decays (BR ~ 96%) close to +1. As the light sbottom mass is 
determined by Mq = 571 GeV, also the ti is light; in fact it is 20 GeV lighter than the 
bi and decays with 96% BR to txi2) thus adding to the signal and the net polarization 
of interest. The 62 and t2 are heavy and play no role for our analysis. The tops stemming 
from ti decays also have a polarization close to +1, so we expect a large positive net 
polarization, diluted only by boost effects. Indeed, we find Vt = 0.75 for this scenario, 
together with rather large asymmetries Ag^ ^ 0.6 and A^^ ^ 0.8. 

• bm-2 is an example of the higgsino scenario with an almost pure RH 61, leading to a top 
polarization close to —1 in 61 -^ txi decays. The ti can be chosen heavy in this case, 
so that the only relevant cross section is 6^6* production. Compared to bm-1, the total 
signal is further reduced by the smaller BR(6i — )■ txi) of only 34%. In the lab frame, we 
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bm-1 


bin-2 


bni-3 


bm-4 


bin-5 


bm-6 


bm-7 


a{pp -^ b^bl) [pb] 


0.137 


0.137 


0.137 


0.137 


0.137 


0.137 


0.137 


a{pp — )■ 62^2) [pb] 


< 10-3 


0.001 


0.082 


0.096 


< 10-3 


0.019 


0.093 


a{pp -^ tit*) [pb] 


0.163 


0.002 


0.095 


0.114 


0.157 


0.033 


0.115 


BR(6i ^txi) 


0.96 


0.34 


0.31 


0.35 


0.72 


0.32 


0.41 


BR(6i ^txs) 


- 


- 


- 


- 


0.16 


- 


- 


BR(62^txr) 


0.30 


0.82 


0.96 


0.96 


0.29 


0.96 


0.58 


BR(62-^tx2) 


<io-3 


0.09 


- 


- 


0.01 


- 


- 


BR(ti^tx?) 


0.51 


0.33 


0.50 


0.50 


0.27 


0.47 


0.01 


BR(ti->tx^) 


0.45 


0.35 


0.47 


0.46 


0.40 


0.42 


0.26 


Vtih^txi) 


0.92 


-0.92 


-0.56 


-0.05 


0.99 


-0.98 


-0.94 


Vt{h^tX2) 


- 


- 


- 


- 


-0.29 


- 


- 


Vt{b2^txi) 


-0.99 


0.99 


0.95 


0.94 


-0.99 


0.99 


-0.97 


Vt{b2^tX2) 


-0.94 


-0.67 


- 


- 


-0.99 


- 


- 


Vtih^txl) 


0.82 


0.40 


0.88 


0.86 


0.99 


0.94 


-0.99 


Vtih^tx'2) 


0.92 


0.50 


0.97 


0.96 


0.99 


0.99 


-0.95 


Vt (total) 


0.75 


-0.73 


0.46 


0.64 


0.92 


0.07 


-0.83 


M 


0.59 


0.14 


0.61 


0.60 


0.80 


0.47 


0.12 


^<Pi 


0.84 


0.57 


0.83 


0.84 


0.92 


0.76 


0.55 



Table 2: Production cross sections at ^/s = 14 TeV, decay branching ratios, top polarizations 
originating from different decay processes in the respective squark rest frame, total polarization 
in the laboratory frame, as well as polar and azimuthal angle asymmetries for the benchmark 
points defined in Table 1. 



have Vt = —0.73 for this scenario, with a small polar angle asymmetry of Ag^ 
a still sizable azimuthal angle asymmetry of A^^ ^ 0.6. 



0.1 but 



bni-3 and bin-4 illustrate the case of intermediate polarization, which is achieved by giving 
the mostly RH 61 a small bi component. Point bin-3 has cos 6^ = 0.1 and Vt (&i -^ tx~[) = 
—0.5, while bm-4 has cos^j^ = 0.15 and Vt (&i — > txi) ~ 0. In contrast to the two previous 
points, both 62 and ti have masses around 700 GeV {i.e. not far from the bi) and thus also 
significantly contribute, in particular because both BR(62 — > txi) and BR(ti — )■ txi2) 
are very large, around 96%. The net polarization in the lab frame turns out to be 0.46 
(0.64) for bm-3 (bm-4). The polar and azimuthal angle asymmetries are sizable, Ag^ ^ 0.6 
and A^j ^ 0.8, for both points. Note however that with roughly 210-250 fb the total rate 
(crxBR's) before cuts for ti + EJp^^^ is smaller than for bm-1 (roughly 290 fb). 

bm-5 is an example for the situation when both charginos are light enough to allow 
bi -^ tXi and 61 — ?> t x^ decays. As before, the LSP is mostly higgsino and the xf ~ Xi 
mass difference is small, such that the Xa ^ Xi ^* is probably missed. The X2 o'^ ^'^^ 
other hand is wino-like and decays to Xi plus an on-shell W . The top polarization in 
the sbottom rest frame is ~ +1 for b\ decays into txi-, while it has a small negative 
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value, Vt = —0.29, for the decays into tx2- As for bm-1, the ti is also light and decays 
with ~ 67% to txi2- The top polarization in these decays, which can likely not be 
distinguished from 61 — )■ t x^, is close to +1. We end up with Vt = +0.92 in the lab frame 
and large asymmetries of A0^ ?^ 0.8 and A^^ ?^ 0.9. 

• bm-6 is a variant of bm-3 with somewhat heavier ti and 62- It is constructed such 
that the net polarization resulting from the combination of Vt{hi -^txi) ~ +1 and 
Vt {ti -^ txi2) ~ "■! is almost zero. Nonetheless the asymmetry parameters are sizable, 

Ae, ^ 0.5 and A^^ ^ 0.8. 

• bm-7 shows the wino-like Xi case, which gives Pj ?^ — 1 in sbottom decays whatever is 
the sbottom mixing angle. Note that the ti decays into txi or tx2 a-lso give Vt ~ —1. 
However, these decays can in principle be distinguished, as m^o ~ 250 GeV and X2 ^ ^Xi 
and Xi — )■ Wxi- To estimate Vt and the leptonic asymmetries, we sum over 61 and 62 
production followed by decays into tx^ . We find Vt = —0.83, Aq^ ^ 0.1 and A^^ ^ 0.5. 

6.3 Polarization dependent kinematic distributions 

In this section we present results for the polarization-dependent kinematic distributions of the 
top decay products introduced in Section 4. As mentioned, these distributions depend on the 
top polarization in the sfermion rest frame, the top boost in the laboratory frame and the 
transverse momentum of the top quark. As a result they are influenced by the mass difference 
between the decaying sfermion and the daughter electroweak-ino as discussed in ref. [22]. Here, 
we concentrate on the case of large mass difference. More concretely, we analyze the lepton 
energy distributions as well as the energy ratios z and u of eq. (37) for bm-2, bm-5 and bm-6, 
corresponding to close to ±1 and zero net polarization, respectively. 

Figure 7 shows the lepton energy Ei (left) and transverse momentum pj, (right) distributions 
for bin-2, bin-5 and bm-6. As can be seen from eq. (33) for a positively polarized top, the decay 
leptons go preferentially in the forward direction and hence, including the boost, the energy and 
transverse momentum of the leptons are larger than in the unpolarized or negatively polarized 
case. For the negatively polarized top, the situation is opposite, the leptons gain less energy 
and the distributions peak at smaller values. This effect is clearly seen in both the Ei and p!p 
distributions in Fig. 7 when comparing the lines for Vt = +0.92, 0.07 and —0.73. 

We conclude that for bi ~ b^ decaying into a higgsino-like chargino (and/or ti ~ t^ decaying 
into top+neutralino [22]), the harder Ei and p-f distributions can be used to extend the reach 
of searches for 3rd generation squarks. On the other hand, for scenarios which lead to small or 
negative top polarization, the opposite conclusion holds. Note, however, that the polarization- 
dependence of the El and p!p distributions makes a general interpretation in terms of SMS [36,37] 
difficult if information on the polarization is used in the analysis. 

The behavior of the Ei and p^ distributions for different polarizations seen in Fig. 7 is a 
reflection of the angular distribution of the lepton w.r.t. the top quark spin direction in the rest 
frame given by eq. (33). Since ki and Kb have opposite signs, the corresponding distributions for 
6-jets will show the opposite behavior. Figure 8 shows the 6-jet energy E^ (left) and transverse 
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Figure 7: Lepton energy Ei (left) and transverse momentum p^ (right) distributions for three 
different polarizations. The dashed blue lines are for bm-2 with Vt = —0.73, the full red lines 
are for bm-5 with Vt = 0.92, and the dotted black lines are for bni-6 with Vt = 0.07. 



Pf = 0.92 
P^ = 0.07 
— - P, = -0.73 




P, = 0.92 
Pi = 0.07 
— - P, = -0.73 




Figure 8: Same as Fig. 7 but for 6-jet energy E'' (left) and transverse momentum p|. (right) 
distributions. 

momentum p^ (right) distributions for bm-2, bm-5 and bm-6. Although the size of the effect is 
smaller then for leptons due to the smaller value of Hb, the &-jet distributions may still provide 
interesting complementary information. In particular, the distributions get harder as the top 
polarization changes from +1 to —1. Thus the loss of reach due to a softened lepton spectrum 
in case of a negatively polarized top might be compensated to some extent by the harder 6-jet 
spectrum. 

In Section 4 we also discussed that the azimuthal angle asymmetry A^p^ and the polar angle 
asymmetry Ag^ may also give us a quantitive measure for the top polarization. To illustrate this 
point, we show in Fig. 9 the azimuthal angle (pi (left) and polar angle 6i (right) distributions 
of the decay leptons for bm-2, bm-5 and bm-6. As expected, the distributions for 0; peak at 
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Figure 9: Distributions of the azimuthal angle (pi (left) and the polar angle 6i (right) of the 
decay lepton. The dashed blue lines are for bni-2 with Vt = —0.73, the full red lines are for 
bni-5 with Vt = 0.92, and the dotted black lines are for bm-6 with Vt = 0.07. 

(pi = and (pi = 27r for all the three cases, but the peaks are higher for a positively polarized 
tops as compared to unpolarized or the negatively polarized ones. A similar situation is seen 
for the distribution of the polar angle 6i: the peaking is again in the direction of the top boost 
and increases when going from negative to positive polarization. The specific values of the 
asymmetries Ag^ and A^^, defined by eqs. (35) and (36), are given in Table 2. The values of 
both the asymmetries are the lowest for negative polarization and increase as the polarization 
goes to 1, thus making them a measure of the polarization. 

The above observables concerning angular distributions of the leptons are independent of 
the anomalous tbW coupling. At LHC-14, however, the tops will be highly boosted and it may 
not be easy to use the angular observables. The boost distribution of the tops for the three 
benchmark points under consideration is shown in Fig. 10. The boost is of course independent 
of the top polarization — the small differences in the boost distributions arise from the different 
stop and sbottom masses for the three benchmark scenarios. The main point is that the tops 
are typically highly boosted. In such a situation, the energy ratios u and z defined in eq. (37) 
may give very useful information. (Recall that these observables can be however affected by 
nonzero values of an anomalous tbW coupling.) A cut on the top boost /3j can enhance the 
dependence on the polarization. 

Figure 11 shows the distributions of the energy ratio u for different cuts on Pf The dis- 
tributions are weighted towards smaller values of u for negatively polarized tops, and towards 
higher values of u for a positive polarization. The cut on the boost enhances the separation 
of +ve and — ve polarization. The analogous distributions of the energy ratio z are shown in 
Fig. 12. As expected from the discussion in Section 4, the behavior of z is opposite to that of u, 
that is positive polarization favours low z values. A cut on Pt again helps to better differentiate 
between different values of polarization. Clearly, asymmetries similar to Aq^ and A^^ can be 
constructed for u and z and may serve as an additional measure of the polarization. 
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Figure 10: The top boost distribution at 14 TeV, for bni-2 with Vt = —0.73 (dashed blue hue), 
bm-5 with Vt = 0.92 (full red hue), and bm-6 with Vt = 0.07 (dotted black hue). 

7 Conclusions 

Given the increasingly strong limits on squark and gluino masses from searches at the LHC, 
together with a SM-like Higgs boson with ruh ~ 125 GeV, "natural SUSY" with light 3rd 
generation squarks, maximal stop mixing, and perhaps light higgsinos emerges as the new 
phenomenological paradigm for SUSY phenomenology to be addressed by the next phase of LHC 
running at 13-14 TeV. It is therefore particularly interesting and timely to develop methods to 
a) enhance the discovery potential for stops and sbottoms and b) determine their properties at 
the LHC. 

In this paper we discussed the polarization of top quarks stemming from sbottom and stop 
decays as a useful tool to this end. In particular we investigated in detail the behavior of the 
top polarization in sbottom decays into charginos, b — ?■ tx~- As in the case of stop decays 
to tx°, this polarization may give clues to the nature of both the sbottom and the chargino. 
Concretly, we pointed out that for a mostly wino-like chargino, the top polarization is always 
fa —1, while for a higgsino-like chargino it can take any value between —1 and +1, depending on 
the sobttom mixing angle. Moreover, we discussed that in realistic setups the relation between 
top polarization and underlying MSSM scenario may not be straightforward because the tops 
may come from several different channels. For example, when the lighter sbottom bi has a large 
LH component, its mass is similar to that of the lighter stop ti (as both are determined by the 
same soft mass parameter), and therefore both stop and sbottom decays may lead to the same 
signature, and only the net polarization resulting from all different decays will be measurable. 

We illustrated the relation between the top polarizations from specific decay channels (in 
the respective stop or sbottom rest frame) and the net polarization in the laboratory by means 
of seven benchmark points. Furthermore, we studied how the top polarization affects its decay 
kinematics and pointed out the strong correlation between the polarization of the top quark 
and the kinematical distributions of the various decay products. We observed how these may 
be used to enhance and evaluate the discovery reach, as well as to construct measures of this 
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Figure 11: Distributions of the energy ratio u, without cut on the boost (top), for (3t > 0.8 (lower 
left) and for Pt > 0.9 (lower right). The dashed blue lines are for bm-2 with Vt = —0.73, the full 
red lines are for bm-5 with Vt = 0.92, and the dotted black lines are for bm-6 with Vt = 0.07. 
The distributions show a strong dependence on the values of polarization considered. 

polarization using the angular observables of the decay lepton. 

In summary, we have shown that top polarization may provide a useful tool for the searches 
for 3rd generation squarks, in particular in the context of natural SUSY with light higgsinos, a 
scenario which is very difficult to resolve at the LHC 
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Figure 12: Distributions of the energy ratio z, without cut on the boost (top), for f3t > 0.8 (lower 
left) and for /3t > 0.9 (lower right). The dashed blue lines are for bm-2 with Vt = —0.73, the 
full red lines are for bm-5 with Vt = 0.92, and the dotted black lines are for bm- with Vt = 0.07. 
The distributions show a strong dependence on the values of polarization considered. 
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